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ABSTRACT 


The unsteady aerodynamic forces and moments of an oscillating airfoil for the 
fixed wing case were determined by Theodorsen along with the development of a lift 
deficiency function. Loewy subsequently developed an analogous lift deficiency 
function for the rotary wing case in which there are an infinite number of layers of 
shed vorticity, or wakes, below the reference airfoil. With the advent of computer 
panel codes that calculate the time histories of the wakes generated by oscillating 
airfoils, a theory is developed for the rotary wing case in which there are a finite 
number of layers of shed vorticity below the reference airfoil. This theory includes a 
lift deficiency function that is completely analogous to Loewy and Theodorsen. 

It has long been recognized that an airfoil oscillating in pure plunge produces a 
propulsive force ("“Katzmayr effect”). Garrick used Theodorsen’s work to develop 
equations for the propulsive force that include the lift deficiency function as a 
parameter. When either Loewy's lift deficiency function or the finite wake lift 
deficiency function is used, the effect of the propulsive force is greatly enhanced with 
the proper phase relationship of the wakes. The finite wake theory along with 
Garrick’s work is used to describe the performance characteristics of Higher Harmonic 
Control. Specifically for the OH-6A, coupled pitch-plunge motion results in a 


propulsive force that significantly reduces the rotor drag force. 
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I. INTRODUCTION 


A. GENERAL 

The present knowledge of rotary wing unsteady aerodynamics is limited. The 
fundamental closed form solutions by Theodorsen [Ref. 1] and Loewy [Ref. 2] provide 
the basis for the theoretical work in this area. While Loewy’s work on wake-induced 
flutter helps explain the phenomenon, it also points out the difficulty to be overcome 
which is the closed form theory shows rapid changes in the lift deficiency function 
with changes in reduced frequency (k), wake spacing (h), and frequency ratio (m). 

In the past, emphasis in the study of unsteady aerodynamics has tended to focus on 
flutter instability and the effects of unsteady aerodynamics on lift generation and 
torsional loads, and not on performance. In this thesis, the emphasis will be on 
performance and the effect of unsteady aerodynamics on the drag of the airfoil. The 
classic work on this subject was performed by Garrick [Ref. 3]. In this work, Garrick 
found that the drag force could be reversed, acting to propel the airfoil forward. This 
drag reversal is generally referred to as either negative drag or propulsive force. 

The primary objective of this thesis is twofold. First, the effect of rotary wing 
unsteady aerodynamics on the drag of an airfoil will be analyzed. It will be shown 
that this effect can be either substantially greater than that shown by Garrick, or 
substantially less depending on the frequency ratio (m). This change in the effect is 


due to layers of shed vorticity (or wakes) below the rotor blades. Second, with the 


recent advent of computer panel code. c: »able of calculating the time histories of the 
wakes generated by oscillating airfoi:-. a :losed form theory is developed for the 
rotary wing aerodynamics in which t. ..re ue a finite number of wakes below the rotor 
blades. It will be shown that this th: ry is completely analogous to Theodorsen and 
Loewy's work, and applicable to Gar. :k: work. 
B. SCOPE 

Chapter II is a short review of th: wc k of Theodorsen [Ref. 1] and Garrick [Ref. 
3]. The material is presented in ari: ine to allow one to understand the key 
assumptions made in each theory w:t!iout being overcome by the details. Chapter III 
contains an analysis of Loewy’s wor: [Re f. 2] and applies it to Garrick. The closed 
form finite wake theory is also presented here, as well as the special case of a single 
wake. Chapter IV is a limited discussion on applications of finite wake rotary wing 
unsteady aerodynamic theory. The first section discusses logical extensions of this 
theory to include forward airspeed [Re!. 4] and compressibility effects [Ref. 5 and 6}. 
A comparison of this theory to the Nava! Postgraduate School (NPS) Unsteady Panel 
Code [Ref. 7 and 8] for the case of pure plunge is also included in the first section. 
The second section presents an explanation for the performance benefits due to Higher 
Harmonic Control (HHC). Appendix A is another application of this theory describing 
the efficient nature of flight of the hummingbird. The discussions in Chapter IV and 
Appendix A are not designed to provide all the answers to the questions arising from 
these applications. They are designed tc show that positive performance effects can be 


found in the application of rotary wing unsteady aerodynamics. 


H. BACKGROUND 


A. THEODORSEN'S LIFT DEFICIENCY FUNCTION 

The aerodynamic forces and moments on an oscillating airfoil for the fixed-wing 
case as determined by Theodorsen [Ref. 1] are based on potential flow theory and the 
Kutta condition. Potential flow theory reduces to Laplace’s equation for the case of 


incompressible subsonic flow. 


2 2 
E (1) 


da oy? 

Since potential flow theory follows the principle of superposition, the potentials can be 
treated as two types: non-circulatory and circulatory. The non-circulatory terms are 
primarily from the effects of the freestream with the airfoil. The circulatory terms 
result from the vorticity generated in the wake of the airfoil. It is in the circulatory 
terms that the lift deficiency function is defined. 

The aerodynamic forces and moments on a thin airfoil, depicted in Figure 1, are 
obtained by integrating the pressure difference across the airfoil determined from the 
generalized Bernoulli equation as follows: 


l 
P=b | p(x)dx (2) 


1 





trailing 
edge 


All distances are non-dimensionalized by the semi-chord (b) 


Figure 1. Parameters of the thin airfoil. 


and 


l 
M, = of p(x)\(x- a) dx , (3) 
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where the generalized Bernoulli equation is written as 


_ n| YOQ , 99 
^ p(x) = 2p (5 - - . (4) 


In the manner suggested by Scanlan and Rosenbaum [Ref. 9], the potential flow 


for the non-circulatory terms can be written as 


Go gh xE-/1-x?J1 zl. (5) 
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The pressure difference across the airfoil for the non-circulatory terms is 
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where the downwash! w(E) = h + b(E-a)a& + va, and simple harmonic motion is 
assumed, or w 7 w,e'?', Integrating the pressure difference in equation (6) to obtain 


the non-circulatory forces and moments yields 


Py, 7 -mpb?(và*h- bad) (7) 


and 
(Mn = ~mp b'|(-- a) vba + D9(_+a?)a-abl (8) 


The potentials of the circulatory terms are determined by integrating the velocity 
potential resulting from a single vortex element with a vortex distribution per unit 
length, U(x,), that is assumed to vary sinusoidally. This yields a circulation potential 


of 


"Only two degrees of freedom (pitch and plunge) are presented here since most rotary 
wing aircraft do not have ailerons incorporated on the rotor blades. 
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The Kutta condition, which requires the velocity at the trailing edge to be finite, can 


be written as 
0 
(22 O (10) 


Ox edge 
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This results in 








w(£) a E dl U(x,) zu dx, . (11) 
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Now defining a lift deficiency function as 
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the Kutta condition becomes 
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Using equation (4) the pressure difference across the airfoil for the circulatory terms: is 
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Integrating the pressure difference in equation (14) to obtain the circulatory forces and 


moments yields 
4 1 y 
mp vb C(X)| va * hs b(7 - à (15) 
and 
] > l ; 
(M). = 2p vb^(a* =) C(&) va * h* aja (16) 


The total pressure difference is obtained by adding the pressure differences for the 
non-circulatory and circulatory terms. After some algebraic manipulation the total 


unsteady pressure distribution can be written as 
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Likewise the total forces and moments are obtained by adding their two parts together 


yielding 


P = -mpb?(và » h- baà) - 2mp vbC(K)| va + he bl, - à] (18) 


and 
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The analytical solution to the lift . efi iency function involves the use of complex 


Bessel functions (or Hankel functio:.s , a: d can be written as 


He? 


CUL) pe ecc (20) 
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where H,” = J, - iY, is the Hankel function of the second kind of order n. The lift 


deficiency function can be separated into its real and imaginary parts as 


C(k) = F(k)+ 1G’ , (21) 





where 
Ho o AA OAE (22) 
C 
and 
Y, + 


AN 
and all Bessel functions are evaluated at the reduced frequency (k). A typical plot of 
F(k) and G(k) as given by Theodorsen [Ftef. 1], is shown in Figure 2. Figure 3 is a 


semi-logarithmic plot of Theodorsen’s lif deficiency function, and the limits of F(k) 
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Figure 3. Semi-logarithmic plot of Theodorsen's lift deficiency function. 


and G(k) as k approaches zero and infinit can be more easily seen. It is an 


interesting note that when E = 0, k = -G = 0.188773655.... 


B. GARRICK'S PROPULSIVE FORCE EQUATION 

One of the earliest investigations on the effect of an oscillating airfoil on the 
resultant drag was conducted by Katzmayr [Ref. 10] in two series of experiments. In 
his second series, Katzmayr measured a 1 egative drag or propulsive force for certain 
velocities and frequencies. This discovery, later to be known as the “Katzmayr 
effect”, helped answer questions or how birds propel themselves through the air. 
Richardson [Ref. 11] used the “Katzmayr effect” to explain the locomotion of fish 
through the water. Garrick [Ref. 3] applied the equations of motion determined by 
Theodorsen [Ref. 1], equations (18) and (19), to a method developed by von Karman 
and Burgers [Ref. 12] to obtain a closed form analytical solution for the propulsive 
force generated by an oscillating airfoil.’ 

Garrick determines the average horizontal force by two methods: energy and 
force. He then compares the results of these methods, and shows they are identical. 
For simplicity only the force method is presented here. First he assumes simple 


harmonic motion for each of the degrees of freedom, resulting in 


@ = a, e 99 (24) 


"Though Garrick developed his equations for three degrees of freedom, only two 
degrees of freedom (pitch and plunge) will be presented here. 
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and 
e E (25) 
where (p, and (p, are phase angles between the initiation of the input and an arbitrary 


reference point. 


The total horizontal force on the airfoil is given by 


P, - wpS?+aP (26) 


where P, is the propulsive force on the airfoil (positive to the left), P is the resultant 
force on the airfoil (positive down), and S is the component of induced downwash in 
the x-direction. The time averaged propulsive force is found by 
= VL d EO 
P. = zl P dt. (27) 
27 Jo 


The x-component of the induced downwash is 


S = gefrati 8 -4à]- bà) . (28) 


In finding the solution to the time averaged propulsive force equation, it is 
arbitrary whether one chooses to employ the real or imaginary parts; the results will be 
the same. Garrick [Ref. 3] chose the imaginary parts, and determined the time 


averaged propulsive force to be 
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It is interesting to note that in the special cases of only one degree of freedom, the 
effects of the arbitrary phase angles between initiation and the reference point 


disappear. The propulsive force equations become for pure plunge (a, = 0) 


X 


P. = apbo! I; (F^ « G?) , (30) 


and for pure pitch (h, = 0) 


B vpbo' bas ([F* « G*]— +(5 -a) ]+5(5 4) (31) 
Fla L) (59 2)2) 


If it is assumed that pitch and plunge are in phase with each other (they reach their 
maximum deflections at the same time), the coupled propulsive force equation 


becomes 


“The result presented here differs from Garrick's equation (34) in reference 3. When 
equation (29) in reference 3 is simplified using the definitions of a, and b,, the result is 
equation (31) above. This result has been confirmed by Garrick's work in reference 10. 


1» 


P, = rp bw? ee + G?) 


1 C 
AAA ARI e» 


+ bah, (> SE lo „S Ąl 





Since propulsive force is analogous to negative drag, it is convenient to define a 


propulsive force coefficient as 


P P 
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The propulsive force coefficients for the special cases become for pure plunge 


(92,7 O) 

Gl (34) 
where h, is the amplitude of deflection non-dimensionalized by the semi-chord, for 
pure pitch (h, =0) 
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and for coupled pitch and plunge 


C, = Tk DE G?)+ai|(FP+ c? y (5 -aP j (5 - a) 2 
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Il} «¢ NALYSIS 


A. LOEWY'S LIFT DEFICIENC* FUNCTION 
The aerodynamic forces and mor. nts on an oscillating airfoil for the rotary-wing 


case are more complex than the fixe< .i: 3 counterpart. Loewy [Ref. 2] determines 


ob 
Qr’ 


different approach to the problem b': -5!v 3g the integral downwash equation which 


these forces and moments for the hos rc. se where k = but uses a slightly 

leads to an equation for the pressure distri dution which is in the same form as 
Theodorsen. Loewy then defines a n:ccif 2d lift deficiency function, and states that 
the integration of the pressure distribution across the airfoil would be identical to 
Theodorsen. This makes the expressicns for forces and moments identical except for a 
modified lift deficiency function. 

The first step that Loewy performs is to set up a system to account for the wakes 
generated by the previous blades in the same revolution and all blades in previous 
revolutions as shown in Figure 4. He uses two indices to account for the vorticity 
shed by a given wake: n indicates the revolution and q indicates the blade whose 
wake it is. The induced velocity or down wash resulting from an element of vorticity 
is obtained from the Biot-Savart Theorem, 


OM O Yy dt 


— n (37) 
2-0) + (1049) 4 
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n = revolution number q = blade number 





Figure 4. Aerodynamic model for a multi-blade rotor system. 


where Y, is the vorticity, Q is the number of rotor blades, and h is the non- 
dimensional wake spacing. Writing the integrals involving the bound vorticity and the 
vorticity in the wake of the reference airfoil separately from the rows of vorticity 


below the plane of the rotor yield 


d UN UE E E ER 
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In essence the problem has been broken down in the same manner as Theodorsen. 


The first integral represents the effects of the freestream on the airfoil (non-circulatory 


I5 


terms). The second integral represents th: velocity created by the vorticity generated 
by the reference wake (circulatory term). The third and fourth integrals represent the 
velocity created by the vorticity generated by previous blades or in previous 
revolutions (circulatory terms). The main difference between Loewy and Theodorsen 
is the terms which account for the vorticity generated by previous blades or in 
previous revolutions. 


Loewy [Ref. 2] shows that the vortici y shed by the qth blade in the nth revolution 


LE ei e At -27mq/Q-27mn)| (39) 
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where I" is the total circulation around the airfoil, y, is the phase angle by which the 
motion of the qth blade leads that of the reference blade, and m is the ratio of 
oscillatory frequency to rotational frequency. Substituting the vorticity expression 


from equation (39) into the integral downwash equation (38) yields 


l caer 
Wa) = Alf Va ar] eT dt 
y x6 


25 X 
Q-1 -. x s AE y — 
e ay e rmalo- V) ND e aun CONS CEN (40) 
q=1 n=0 -00 (x-Ey $ (nQ* gy 4? 


nz] Le (x-EY + n? Q?p? 


The last two integrals in equation (40) have the form 


= (x- Ode 2 ime tira (41) 
-0% (x- EY + A? 
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Substituting equation (41) into equation (40) and noting that the summations over n 


are convergent geometric series give 





i a 
d yY,(5) dG ie ee eK dé 5 -ikx , (42) 


where 
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The function W(k,h,m) may be thought of as a weighting function for the vorticity 
shed by previous blades or in previous revolutions. 

The form of the downwash equation in equation (42) can be solved by applying 
Sohngen’s inversion formula [Ref. 14], which shows that the solution to an equation in 


the form 


Abd (44) 


g(x) = = 
= 





Ax) = E = af H As | (45) 
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Satisfying the condition f(1) = finite is the same as employing the Kutta condition. 


The bound vorticity becomes 
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Evaluating the circulation over the entire airfoil yields 
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where the Hankel and Bessel functions are evaluated at reduced frequency (X). 
Since the airfoil can be thought of as a vortex sheet, the generalized Bernoulli 


equation (4) becomes 
Q X 
apo) - -p|vy,cn * 5 [ (E DE). (48) 
=l 
If simple harmonic motion is assumed, 
£ 12 = 10 Y; 2 (43) 
Ot 
and equation (48) becomes 


Ap(x) - -p|vy,Gsn) + Io of 15,0 dE (50) 


Substituting the bound vorticity equation (46) and the airfoil circulation equation (47) 


into equation (50) yield 
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The pressure distribution in equation (51) has the same form as Theodorsen's in 


equation (17) if the factor multiplying the first integral is written as 
i 2 C's hm) , (52) 
T 

where C(k,h,m) is Loewy's lift deficiency function. Solving for C'(k,h,m) yields 


HO «2J, Wh 
C'(khm) = ee : (53) 
Hu T * 217, * id W(k, hm) 


where the Hankel and Bessel functions are evaluated at reduced frequency (k). Since 
C'(K, h, m) is not a function of location along the airfoil, the integration of the pressure 
distribution in equation (51) across the airfoil will yield equations of motion identical 
to Theodorsen’s in equations (18) and (19) except C(k) will be replaced by C (k, h, m). 
It can be seen that as W(k,h,m) approaches zero, C’(k,h,m) = <(k), and this condition 
corresponds to an infinite wake spacing (h-o). 

Since the wake weighting function is periodic in nature, Loewy [Ref. 2] shows that 


the wake weighting function for a multi-blade rotor can be expressed by a single-blade 
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rotor with modified values of h and » th - yield the same value of W. For a single- 


blade rotor, the wake weighting func’ »n ‘ecomes 


Wu CE (54) 


»Kh ¿12m E 


where A and m are the modified val.. > oœ: wake spacing and frequency ratio 
respectively. For the case of the sir¢ 2-b ade rotor the wake weighting function can 


be separated into its real and imagi! : p. rts, 


WA hiv) = a+ iB , (55) 
where 
: cos 2m - e ^^ 
E E a (S6) 
e*^ . 2cos2mm » e * 
and 


sin 277 12 
B.——— ___. (57) 
e*” - 2cos 21m + e 


Separating Loewy's lift deficiency function into its real and imaginary parts similar to 


Theodorsen yields* 


Ck bm) = F'+iG' (58) 


where 


“This is equivalent to Loewy’s equations in reference 2 when the substitutions 
a = 1+ 24 and B = -26 are used. 
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. 402244 - 24) B (59) 


F' 
A? + B? 
gi - M2 ABA +A +24)B (60) 
A? + B? 
and 
A = J(1+26) + Y, -2JB (61) 
B= -Y¥,+2/8+4(1+24) . (62) 


Plots of Loewy's lift deficiency function are shown in Figures 5 through 8. It is easily 
seen that there are rapid changes in the lift deficiency function as reduced frequency 
(k), wake spacing (h), and frequency (m) vary. 

Since Loewy's lift deficiency function is completely analogous to Theodorsen's lift 
deficiency function, the propulsive force determined by Garrick [Ref. 3] and the 
propulsive force coefficient can be calculated for the rotary-wing case by substituting 


F’ and G’ for F and G. For pure plunge (a, = 0) 
Cp = Th, (F’)? + (G’)?) (63) 
for pure pitch (h, = 0) 


G, Rd -ailh a) 


1 1 ] E 
a A. 


(64) 
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Figure 5. Loewy's lift deficiency function as a function of wake spacing for 
m = 0. 
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Figure 6. Loewy’s lift deficiency function as a function of wake spacing for 
m = 0.25. 
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Figure 7. Loewy’s lift deficiency function as a funciton of wake spacing for 
m = 0.5. 
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Figure 8. Loewy’s lift dieficiency function as a function of wake spacing for 
m = 0.75. 
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and for coupled pitch and plunge 


C, = 7k D (P ^ 





«as (P? «o L (Ža) y 
1,1 ] ] ] o 
e 

+00, (5-3) y? (G')?)+ (5 — eJ] 


Plots of propulsive force coefficient « «3 inction of frequency ratio for k = 0.1234 are 
shown in Figure 9 through 11. It ci: be : een that Cp, is reasonably constant for 0.15 
< m < 0.8 and greater that the “Katz yr affect.” (In this case h = 100 is sufficient to 
represent h-). For m < 0.15 and m > 0.8, Cp, begins to decrease to the point where 


it eventually becomes less than the “h:.izmayr effect.” 


B. FINITE WAKE LIFT DEFICIENCY FUNCTION 

1. General 

Loewy's theory can be extenc«d to the case of a finite number of wakes 

below the rotor blade by modifying the wake weighting function. This modified wake 
weighting function is completely analczovs with Loewy's which makes the definition 
of the lift deficiency function the same, but with & and f modified by the new wake 
weighting function. Since the lift deficiei cy function is of the same form, the 
aerodynamic forces and moments wil! be the same, but with a modified lift deficiency 
function. Loewy has already shown that ‘he rotary-wing problem reduces to the case 


of a single rotor. The modified wake we: zhting function will also be developed for 
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Figure 9. Propulsive force coefficient in pure plunge as a function of wake spacing 
with an infinite number of wakes (Loewy’s lift deficiency function). 
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Figure 10. Propulsive force coefficient in pure pitch as a function of wake spacing 
with an infinite number of wakes (Loewy’s lift deficiency function). 
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Figure 11. Propulsive force coefficient in coupled pitch-plunge as a function of wake 
spacing with an infinite number of wakes (Loewy’s lift deficiency function). 


the case of the multi-blade rotor, and then reduced to the case of the single-blade .rotor 
with modified wake spacing and frequency ratio. 

The development of a finite wake lift deficiency function is the same as 
Loewy's described earlier. Changing the integral downwash equation (40) to a finite 


number (N) of revolutions yields 


w(x) = amr Ya a | pened 
^ x-t 
- HAT me VOX Linn Cranage (66) 
n p AO 


E 2) a 2 e ik (x - Eide 
ENGEL OS 
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where the summations from n equal zero to infinity have been replaced by finite 
summations from n equal zero to N. Using equation (41) to solve the last two 
integrals in equation (66) results in 
d E 
i al ee ae E ar] es qe 
2T me c 


B i N 
z m&T ei pue aa e mne -KaQ- 9h (67) 


g=1 n-0 


N 
+ Ye inane -nQkb | 


n=1 


Since this is in the same form as eguation (42), the terms enclosed by the parentheses 
can be defined as the finite wake weighting function, or 
W, (Kh) 2 Y ¿rado ey y e 7n g-kaQ-9)5 , y e i2rma ¿-nQkb (68) 
g n-0 n-1 
It is easily seen that this modified wake weighting function is periodic once the 
relationship y, is known, and a multi-blade rotor can be reduced to a single-blade 
rotor with modified wake spacing and frequency ratio. For the case of the single- 


blade rotor (Q = 1), 


N 
WAKA, m) = Y emma ¿-nkb (69) 
n=l 
It can be shown that as N-~, the finite wake weighting function in equation (69) 
becomes a convergent geometric series equivalent to Loewy's wake weighting function 


in equation (54) provided m + O and h + 0. A singularity in the solution exists when 


m = O and A = O, and the case of A = 0 has no real physical significance. 
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The solution to the remainder of the finite wake problem is identical to 
Loewy, but using the modified finite wake lift deficiency function. To keep the finite 
wake lift deficiency function separate from Loewy's lift deficiency function, the former 


is defined as follows: 


H? +27 W 
C"(kh,m) = ——— 1 -——— , (70) 
FR” + iFy” « 2U, * i4] Wy 


Separating the finite wake weighting function in equation (69) into real and imaginary 


parts, 
Wy = &y* By, (71) 
yields 
N 
EOS Zn (72) 
n-i 
and 
N 
e S (73) 
n=} 


The finite wake lift deficiency function becomes 


E (74) 


where 
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F’ J (I an OAS (75) 
An+ By 
Ge "(Y 224 B AN AA +28 By (76) 
AU CBS 
and 
AS 20 Y 229] 5. (77) 
By = =Y +21B,* 4 (1+26y . (78) 
2. Single Wake Case 
For the case of a single wake, the wake weighting function reduces to 
W = e armo Kb) (79) 
where 
à, - e "^ cos2rm (80) 
and 
(81) 


B, = -esin 2am . 
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Plots of the single wake lift deficienc: f ction are shown in Figures 12 through 15. 
The propulsive force coeffic at an be calculated by replacing F and G in 


equations (34), (35), and (36) with F inc G'. For pure plunge (a, = 0) 


G - TE (GC) (82) 


for pure pitch (h, = 0) 


C, - m Ea (FP (G) la (a) 





7 1 1 (G') t 
xn We ee d, D | , 
and for coupled pitch and plunge 
Cp, = ze Ac Qnae 
: l 2 
«a? (um y «cay Za) 
Ww ui £ 





Plots of the propulsive force coefficient is a function of frequency ratio for k = 0.1234 
are shown in Figures 16 through 18. It is interesting to note that for 0.24 < m < 0.7 
the propulsive force is greater than tlic ")xatzmayr effect" with the maximum occurring 


just prior to m = 0.5. Comparing Figure: 9, 10, and 11 with Figures 16, 17, and 18 
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Figure 12. Single wake lift deficiency function as a function of wake spacing for 
m - Q. 
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Figure 13. Single wake lift deficiency function as a function of wake spacing for 
m = 0.25. 
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Figure 14. Single wake lift deficiency function as a function of wake spacing for 
m = 0.5. 
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Figure 15. Single wake lift deficiency function as a function of wake spacing for 
m - 0.75. 
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Figure 16. Propulsive force coefficient in pure plunge as a function of wake spacing 


with a single wake. 
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Figure 17. Propulsive force coefficient in pure pitch as a function of wake spacing 


with a single wake. 
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Figure 18. Propulsive force coeffici:1:: i: coupled pitch-plunge as a function of wake 
spacing with a single wake. 


respectively, the maximum propulsive force (at approximately m = 0.5) for a single 
wake is always greater than the infinite number of wakes. The reason for this increase 
is that the absolute value of C and C' «e never greater than 1.0, but C as m 
approaches 0.5 exceeds 1.0 for certain wake spacings. 

When m = 0.5, this corresponds to the wakes being 180° out of phase as 
shown in Figure 19, and the vortex from the wake of the reference blade in the current 
revolution is directly above an opposiic!y spinning vortex from the single wake of the 
previous blade or revolution. This efficie it use of the vortices in the wake causes C' 
to exceed 1.0. As the number of wakes increases the vortices from lower wakes are 
not always 180° out of phase with the reference blade. When m = 0.5, each wake is 


180° out of phase with the wake above ad below it. Thus, each wake is in phase 
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Figure 19. Vortex interaction when wakes are 180° out of phase (m = 0.5). 


with the wakes +2 wake spacings from it. Since the finite wake weighting function 
contains the term e*””, the value of each term in the wake weighting function will 
alternate in sign. However, since the wake spacing is greater for each subsequent 
wake, the absolute value of each term is not as large as the first term, and thus, the 
effect of the in-phase wakes are not as strong. As the wake spacing increases this 
effect diminishes. 
3. Multiple Wake Case 

When the number of wakes increases, the value of frequency ratio at which 

the maximum propulsive force coefficient occurs changes as shown in Figures 20 


through 22 for h = 2.0. This change is analogous to the estimation of a square wave 
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by a finite number of terms in a Fourier series, where the propulsive force coefficient 
determined using Loewy's wake weighting function is considered the square wave. As 
Figures 16 through 18 and 20 through 22 show, an increase in either wake spacing (A) 
or the number of wakes (N) will decrease the maximum propulsive force that can be 


obtained. The limit of this decrease will be Loewy's case for a hover. 
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Figure 20. Propulsive force coefficient in pure plunge as a function of the number of 
wakes at h = 2.0. 
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Figure 21. Propulsive force coefficient in pure pitch as a function of the number of 
wakes at h = 2.0. 
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Figure 22. Propulsive force coefficient in coupled pitch-plunge as a function of the 
number of wakes at h = 2.0. 
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IV. APPLICATIONS 


A. EXTENSIONS AND COMPARISON OF THEORY 
1. Forward Airspeed 
Extension of the finite wake theory can easily be made to include forward 
airspeed. Shipman and Wood [Ref. 4] expand on Loewy’s work incorporating forward 
airspeed via the advance ratio (p) into the rotary wing unsteady aerodynamic problem. 


The definition of reduced frequency becomes 


ob 
2 8 
CE SS 


The Biot-Savart Theorem is used with wakes extending out to infinity. Since 
freestream velocity at any point on the rotor blade is a function of both blade section 
radius and azimuth position, the shed vorticity will experience a build-up and decay as 
azimuth varies from 0° to 360°. Shipman and Wood developed a decay function to 
account for the variation in vorticity and incorporate this function and the advance 
ratio into the integral downwash equation. The pressure distribution is found by 
applying Sóhngen's inversion formula [Ref. 14] to the downwash equation. The 
pressure distribution is shown to be in the same form as Theodorsen and Loewy, and a 
modified lift deficiency function is defined to account for the differences. Since this 
modified lift deficiency function includes Loewy's wake weighting function as one of 


its terms, the finite wake weighting function could be directly substituted into Shipman 
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and Wood's lift deficiency function. If the segment of the wake weighting function 
due to the vorticity decay (AW) is modified to account for a finite number of wakes 
instead of an infinite number, a new finite wake forward airspeed lift deficiency 
function can be defined, and the remainder of the problem is the same as Shipman and 
Wood. 

2. Compressibility Effects 

Adding the effects of compressibility to the rotary wing unsteady aerodynamic 
problem was first accomplished by Jones and Rao [Ref. 5]. Hammond [Ref. 6] 
provided an alternate approach to the problem shortly after Jones and Rao, but only 
the method by Jones and Rao will be presented here. 

Jones and Rao expand on Loewy's work by taking the known solution to the 
fixed wing two-dimensional unsteady airfoil theory and modifying it to include the 
effect of an infinite number of wakes below the reference airfoil (or rotor disk). 
Euler's equation is used with a coordinate change to determine the downwash equation. 
An interesting result is found, and that is for a given circulation, the downwash due to 
an infinite system of wakes is the same for compressible and incompressible flow. In 
other words, the circulation generated by layers of shed vorticity beneath an airfoil is 
independent of Mach number. This fact will be true whether there is an infinite 
number of wakes or a finite number of wakes. Thus, the finite wake theory can be 


applied here, and the remainder of the problem is the same as Jones and Rao. 
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3. NPS Unsteady Panel Code 
The NPS Unsteady Panel C> 2 v 4s written by Pang [Ref. 7] to solve the 
potential flow for two airfoils execut. 3 u.isteady motion in an inviscid, 
incompressible flow. In this code, t!. tw > airfoil surfaces are approximated by a large 
number of surface elements, and a li: or 1 source distribution and vorticity are placed 
on each element. The source streng ! ve’ es from element to element , while the 
vortex strength is the same for all ese: :er 5, The singularity strengths are determined 
from the flow tangency condition or. < >t}. airfoil surfaces and the two Kutta 
conditions. 
Extensive comparisons of th:: co te with Theodorsen was accomplished by 

Riester [Ref. 17], and the code (using an NACA 0007 airfoil) shows very good 
agreement with flat plate theory. Comparison of this code with Loewy could not be 
accomplished easily since it would rec uir? a very large number of airfoils to generate 
the infinite number of wakes. Using t! > finite wake theory, the code could easily be 
compared to the case of a single wak:. A comparison for pure plunge [Ref. 8] with h 
= 2.0, k = 0.1234, and h, = 0.14 are presented in Figure 23. The results show that the 
NPS Panel Code has very good agreement with the finite wake theory. 
B. HIGHER HARMONIC CONTRO! 

Higher Harmonic Control (HHC) i: «n active vibration control system for 
helicopters. The concept of HHC is to a-ter the aerodynamic loads on the rotor blades 


such that the blade response is reduced, vhich in turn reduces the vibratory forces and 


44 


Single Wake Theory 


EN 
| 
| 
! 


Fatzmüyt Effect ' 


1 


| 
! 
| 


! | 


b = Panel Code (2 airfoils 
k = Panki Code (1 airfoil) 

MEAS ioe 
0.3 0.4 0.5 0.6 


Frequency Ratio (m) 





Figure 23. Comparison of propulsive force coefficient in pure plunge (h = 2.0, k = 
0.1234, and h, = 0.14). 


moments acting at the hub that cause vibration. HHC flight test results (Ref. 15], 
under a joint NASA-Army sponsored program using a modified OH-6A, show a 
reduction of airframe vibrations up to 90 percent. Basically, HHC is an electronic, 
computer-controlled active vibration suppression system which senses and cancels 
vibrations in a helicopter by Q per revolution feathering or pitch motion of the rotor 
blades, where Q is the number of rotor blades? On the OH-6A, higher harmonic 
blade pitch control was achieved by superimposing 4/rev (or 4P) swashplate motion on 
basic cyclic and collective control inputs. The aircraft was flown from hover to 100 


knots, but the remainder of the discussion will focus on the hover performance. 


*Normally for helicopters, N is used to represent the number of blades, but to 
remain consistent with Chapter IH Q will be used. 
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In addition to reducing vibration levels HHC on the OH-6A showed the potential 
for decreased helicopter power requirements. Figures 24 and 25 show hover power 
required for three different schedules of HHC open loop excitation. What is of 
interest here is that the amount of power reduction achieved (on the order of 10% or 
20 horsepower) appears to be relatively independent of the type of excitation applied. 
This power reduction is important, especiclly in hover where power required is large, 
because of the benefit of lowering fuel consumption, which can be translated into a 
larger payload or increased mission time. 

The physics behind HHC can be explained by unsteady aerodynamics. The input 
by HHC to the rotor blade is a pitch oscillation, +s degree on the OH-6A. From 
Figure 10 (infinite wakes) it can be seen that the propulsive force coefficient is 
negative (drag) for the wake spacing of the OH-6A (h = 2.0) regardless of the phasing 
of the wakes (m). Thus, it might be concluded that pitch oscillations alone cannot 
account for the substantial reduction of power recorded by the flight tests. If it is 
assumed that viscous interaction of vortices reduce the problem to a finite wake 
solution as shown in Figures 17 and 21, the propulsive force obtained is still not 
sufficient to account for the power reduction even if the wakes were optimally phased. 

As stated earlier, HHC provided a 4P pitch input in both collective and cyclic were 
applied at the blade root. The collective input results in a 4P oscillation along the 
blade, while the cyclic input results in 3P and 5P oscillations along the blade. From 
the Southwell plot of the OH-6A [Ref. 16] shown in Figure 26, it can be seen that the 


second flapwise and third flapwise natural frequencies occur near 3P and 5P 
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Figure 24. Effect of HHC on main rotor shaft power in hover for selected set of three 
open loop inputs. 


respectively. During the course of a revolution of a blade about the hub, the blade 
will naturally experience flapping motion due to changes in angle of attack, or pitch 


changes, as the blade proceeds around. Since 3P and 5P flapping frequencies are 
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Figure 25. Effect of HHC on main rotor shaft power in hover for selected set of two 
open loop inputs. 


near natural frequencies, small 3P and SP pitch inputs will result in much greater 
flapping outputs. Thus, the unsteady aerodynamics of HHC for the OH-6A is coupled 


pitch-plunge oscillations. Figures 11, 18, and 22 show that coupled pitch-plunge 
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Figure 26. Southwell plot of the OH-6A rotor blade. 
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oscillations always result in a propulsi e orce regardless of the number of wakes 
considered or phasing of those wake. W aat is of greater importance is that if the 
wakes are optimally phased together, .e , ropulsive force is an order of magnitude 


higher than pure plunge or pure pitc! T} > bottom line for the OH-6A is that the 


tv 
PA 


inherent natural frequencies of the b: long with 4P cyclic inputs result in coupled 


pitch-plunge oscillations which yielc . pi oulsive force large enough to account for a 


T3 


substantial power reduction for the 1« 20; er. 
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V. CONCLUSION 


The unsteady aerodynamics of rotary wings is inherently more complicated 
because of the layers of shed vorticity, or wakes, beneath the rotor disk that interact 
with the rotor blades. The Biot-Savart Theorem is used to account for the effects of 
these wakes, and a modified lift deficiency function is developed for use with 
Theodorsen's equations of motion [Ref. 1]. Loewy [Ref. 2] develops a lift deficiency 
function for the case of an infinite number of wakes, while this thesis develops a lift 
deficiency function for any finite number of wakes. This thesis also presents the 
special case of a single wake, and shows that for small wake spacings there are certain 
reduced frequencies where the lift is actually enhanced. In other words, the lift 
deficiency function becomes a lift efficiency function. 

The equations describing the propulsive force generated by an oscillating airfoil 
are developed by Garrick [Ref. 3]. These equations are derived from Theodorsen’s 
equations of motion, and included the lift deficiency function as a parameter. Since 
the rotary wing lift deficiency functions (Loewy’s and the finite wake) are just 
modifications of Theodorsen's work, Garrick's work can be applied to rotary wing 
ur.steady aerodynamics. When either Loewy’s lift deficiency function or the finite 
wake lift deficiency function is used, the propulsive force can be greatly enhanced 
with the proper phase relationship of the wakes. The use of coupled pitch and plunge 


can result in a propulsive force which is an order of magnitude higher than either pure 
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plunge or pure pitch. This coupling of plunge and pitch is used to explain the 
performance benefits reported in the flight test data for the OH-6A equipped with 
Higher Harmonic Control. Coupled pitch-plunge oscillations with proper phase 
relationships of the wakes yield a propulsive force large enough to account for a 


substantial power reduction for the helicopter. 
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APPENDIX A 


I. FLIGHT OF THE HUMMINGBIRD 

For thousands of years, man has been intrigued by the flight of birds. Early Greek 
mythology includes a tale of Daedalus and Icarus who made wings held together by 
wax, and escaped their imprisonment on the island of Crete by flying through the air 
like birds. Early attempts by man to fly were usually designed around imitations of 
birds. Most of these early attempts met disastrous results because of a lack of clear 
understanding of the aerodynamics of ornithological flight, and that is the wing on a 
bird is a coupled lift and propulsion device. In fact, the Wright brothers first 
successful flight demonstrated man's limited knowledge of aerodynamics by 
decoupling lift and propulsion in order to achieve flight. To this day lift and 
propulsion remain decoupled. 

The coupling of lift and propulsion in birds can be explained by unsteady 
aerodynamics. As a bird flaps its wing down, the trailing edge also begins to rise 
[Ref. 18]. Thus, the motion of the wing is a coupled pitch-plunge motion. It was 
shown earlier that coupled pitch-plunge motion yields a propulsive force that is an 
order of magnitude higher than pure plunge. Since most birds require some forward 
airspeed before being capable of achieving flight and the rate of flapping is between 1 
to 5 beats per second, the analogy to unsteady aerodynamics can be described by 


Garrick’s equation for coupled pitch-plunge using Theodorsen’ lift deficiency function 
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(equation (36)). Birds take advantage of the extra propulsive benefit of coupled pitch- 
plunge to make flight appear effortlessly. 

A few birds, most notably the hummingbirds, are capable of hovering. In hovering 
flight the hummingbird’s body is angled upward at about 45°, so that the plane of 
wing flapping is approximately horizontal. As the wing moves up (backward) and 
down (forward), the angle of attack, or pitch, changes making a figure eight pattern 
that is on its side. In essence, the hummingbirds wing motion is similar to the motion 
of the hands of a swimmer when he is treading water. Once again the hummingbird 
takes advantage of a coupled pitch-plunge motion to achieve a higher propulsive force. 

What is unique about the hummingbird is the rate of flapping, or the oscillation 
frequency; it is on the order of 40 to 80 beats per second. Welty (Ref. 18] states that 
the hovering flight of the hummingbird is grossly inefficient when compared to 
helicopter rotor blades, but this is overcome by an increase in apparent velocity 
created by trailing currents from the previous reverse stroke. This is not exactly 
correct. The hummingbird’s wing motion is similar to that of helicopter rotor blades. 
Both are cyclic motions of changes in pitch (angle of attack) and plunge (flapping). 
The hummingbird is just as efficient as a helicopter - both have a large power required 
in hover. The apparent increase in velocity can be explained by trailing vortices. 
Since the hummingbird’s rate of flapping is large in a hover, the wake created by a 
previous stroke will not dissipate and be carried downstream. Each stroke of the wing 
always has layers of shed vorticity directly below it. If these wakes are optimally 


phased depending on the number of wakes, the benefit will be a propulsive force 
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greater than that of non-hovering birds. The hummingbird must use wakes efficiently, 
or it will not hover. Thus, the analogy of the hummingbird to unsteady aerodynamics 
can be described by Garrick’s equation for coupled pitch-plunge using the finite wake 


lift deficiency function, or in the case of a pure hover Loewy's lift deficiency function. 
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